Viscoplasticity and large-scale chain relaxation in glassy-polymeric strain hardening 
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A simple theory for glassy polymeric mechanical response that accounts for large scale chain 
relaxation is presented. It captures the crossover from perfect-plastic response to Gaussian strain 
hardening as the degree of polymerization N increases, without invoking entanglements. By relating 
hardening to interactions on the scale of monomers and chain segments, we correctly predict its 
magnitude. Strain activated relaxation arising from the need to maintain constant chain contour 
length reduces the characteristic relaxation time by a factor ~ eN during active deformation at 
strain rate e. This prediction is consistent with results from recent experiments and simulations, 
and we suggest how it may be further tested experimentally. 

PACS numbers: 61.41.+e,62.20.F-,81.40.Lm,83.10.Rs 



I. INTRODUCTION 

Developing a microscopic, analytic theory of glassy 
polymeric mechanical response has been a longstanding 
challenge. Plasticity in amorphous materials is almost al- 
ways viscoplasticity, i. e. plasticity with rate dependence. 
Many recent studies have focused on plasticity in metallic 
or colloidal glasses. Relative to these systems, polymer 
glasses possess a wider range of characteristic length and 
time scales, because of the connectivity, uncrossability, 
and random- walk-like structure of the constituent chains. 
These alter the mechanical properties significantly [l[ ; for 
example, a uniquely polymeric feature of plastic response 
is massive strain hardening beyond yield. 

The relationships between polymeric relaxation times 
on different spatial scales are fairly well understood for 
melts [31, but much less so for glasses. Recent experi- 
ments [IHH , simulations 0-0] , and theories @, Q have 
all shown that local (segment-level) relaxation times in 
polymer glasses decrease dramatically under active defor- 
mation (especially at yield) and increase when deforma- 
tion is ceased. Analysis of this phenomena has focused 
on stress-assisted thermal activation of the local relax- 
ation processes, but it is likely that other structural relax- 
ation processes at larger scales or of different (c. g. strain 
activated) character are also important in determining 
the mechanical response. Improved understanding of the 
concomitant scale-dependent relaxation is necessary to 
better understand polymeric plasticity and material fail- 
ure. However, theoretical prediction of large-scale relax- 
ation in deformed polymer glasses is still in its infancy; 
most treatments evaluate mechanical response "neglect- 
ing the effect of the ongoing structural relaxation during 
the experiment" [Icl |. 

This paper is an attempt to improve on this state of 
affairs. We develop a theory that treats large-scale relax- 
ation of uncrosslinked chains during active deformation. 
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Stress in the postyield regime is assumed to arise from 
the local plastic rearrangements similar to those which 
control plastic flow; as strain increases, these increase 
in rate with the volume over which they are correlated. 
Polymeric strain hardening is thus cast as plastic flow in 
a medium where the effective flow stress increases with 
large scale chain orientation. Relaxation of chain orien- 
tation is treated here as inherently strain activated and 
coherent, i. e. cooperative along the chain backbone. 

Our theory predicts a continuous crossover from per- 
fect plasticity to "Gaussian" (Neohookean) strain 
hardening as the degree of polymerization N increases. 
The latter form is predicted when chains deform affinely 
on large scales, and corresponds to the limit where the 
system is deformed faster than chains can relax. A 
key difference from most previous theories is that in- 
stead of invoking entanglements, we relate the timescale 
r for large scale chain relaxation to the segmental re- 
laxation time r a . This is consistent with: (i) the pic- 
ture that stress arises predominantly from local plastic- 
ity [lH, (ii) recent dielectric spectroscopy experiments 
indicating connections between relaxations on small and 
large scales [Hj], and (iii) recent NMR experiments [l4| 
that have found the effective "constraint" density for de- 
formed glasses is much larger than the entanglement den- 
sity measured in the melt. By relating strain hardening 
to interactions on the scale of monomers and segments, 
we make a novel prediction of its magnitude. We test 
our predictions using coarse-grained molecular dynamics 
simulations of polymer glasses, and in all cases find (at 
least) semiquantitative agreement. 

The rest of this paper is organized as follows. In Sec- 
tion|II]we motivate and develop the theory for mechanical 
response, make predictions that illustrate the effect of co- 
herent relaxation in constant-strain-rate deformation and 
constant-strain relaxation experiments, and test these us- 
ing simulations. Finally, in Section lnTl we summarize our 
results, place our work in the context of recent theories 
and experiments, discuss how the model could be more 
quantitatively tested experimentally, and conclude. 
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II. THEORY AND SIMULATIONS 
A. Background 

Consider a bulk polymer sample deformed to a macro- 
scopic stretch A. Classical rubber elasticity relates the 
decrease in entropy density to g(X) — | {X 2 X + X y + X 2 ) [H, 
[l6j . Phenomenological 'Neohookean' theories assume a 
strain energy density of the same form. Both approaches 
give an associated (true) stress a oc dg{X)/dh\(X) = 
g(X). All results in this paper are presented in terms 
of true stresses and strains. Stress-strain curves in well- 
entangled polymer glasses are often fit [HI by a(X) = 
Co + Gng(X), where ao is comparable to the plastic flow 
stress <j fi OWl and Gr is the strain hardening modulus. 
Because of this, strain hardening has traditionally been 
associated with the change in entropy of an affinely 
deformed entangled network, with Gr assumed to be pro- 
portional to the entanglement density p e . Forms for <?(A) 
and g(X) for the most commonly imposed deformation 
modes are given in Table HJ Figure Ha) depicts the "uni- 
axial" case. 

There are, however, many problems with the entropic 
description [18l - [20l |. One is that chains in uncrosslinked 
glasses will not in general deform affinely at large scales 
comparable to the radius of gyration. Rather, they will 
possess a chain-level stretch X e ff (Figure [Ha)) which de- 
scribes the deformation of chains on large scales; for ex- 
ample, the zz-component is (Rz/Rl), where R z is the 
z-component of the rms end-to-end distance i? c , and RP Z 
is its value in the undeformed glass. The deformation 
of well-entangled chains is consistent with an affine de- 
formation, X e ff = A, while for unentangled systems the 
deformation is subaffine (2lj . Hereon we drop tensor no- 
tation, e. g. A — > A, but all quantities remain tensorial. 

A key insight is that the evolution of stress in poly- 
mer glasses is controlled by the stretch X e f / of chains on 
scales comparable to their radius of gyration, and only 
indirectly by A. Well below the glas s transition temper- 

mM by 



ature T g , stress is well described 



a(X) =a + G R g(X eff ), 



(1) 



where G R is the value of Gr in the long-chain limit [25| ■ 

Equation[T]shows that predicting A e / / is a critical com- 
ponent of the correct theory for the mechanics of un- 
crosslinked polymer glasses. However, to our knowledge, 
no simple microscopic theory that predicts the functional 
form of X e ff in glasses has been published. Most vis- 
coelastic and viscoplastic constitutive models that de- 
scribe strain hardening, e. g. Refs. [2^, [23] , decompose 
A into rubber-elastic and plastic parts, or use other in- 
ternal state variables, but do not explicitly account for 
A_7^ X e ff or the A^-dependence of nonaffine relaxation 
|2q . In this paper we do so; A e // is treated as a meso- 
scopic (chain-level) internal state variable [28j . 
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FIG. 1: (Color online) Schematic of our model, (a) A is the 
macroscopic stretch (observed on the scale of the experimen- 
tal sample if deformation is homogeneous), while X e ff is the 
large-scale chain stretch. Constant-volume 'uniaxial' defor- 
mation is assumed, (b) The spring-dashpot model for the 
nonaffine chain response is described by Eq. [2] 



TABLE I: Functional forms for strain hardening assuming 
affine constant-volume deformation by a stretch A. 



Def. mode 
Uniaxial 
Plane Strain 
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9(X) 
A 2 — i 



5(A) 
i(A 2 + 1 + ^) 



Maxwell-like model for A 
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We now develop a predictive theory for A e //. The 
problem is most naturally formulated in terms of true 
strains e e ff — m (^e//) and e = ln(A). We postulate a 
simple Maxwell-like model for the relaxation of e e f / . The 
governing equation for the model shown in Figure Wfo) is 



eff/r- 



(2) 
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Equation [2] is a standard "fading memory" form imply- 
ing chains "forget" their large-scale orientation at a rate 
t -1 . In other words, r is the time scale over which e e // 
will relax towards its "equilibrium" value e e ff = (we 
assume that chains are not crosslinked) . 

In this formulation e e j f corresponds to the strain in the 
"spring" in Fig. [ljb). However, the stress we will asso- 
ciate with increasing |e e //| is viscoelastoplastic. Here e e ff 
is an (in principle) micro-reversible strain corresponding 
to chain orientation; it is not an elastic strain in the 
macroscopic sense of shape recoverability of a bulk sam- 
ple [34j . e p i = e — e e ff corresponds to the "dashpot" 
strain used in many constitutive models; it is plastic in 
the sense of being both micro- and macro-irreversible. 

Maxwell-like models have been used to describe poly- 
mer viscoelasticity for more than half a century, and com- 
plicated ladder models were developed (e. g. Ref. (35j ) 
because of the inadequacy of earlier single rate models. 
However, we will provide evidence below that the correct 
choice of meso- variable (i. e. A e // or e e //) restores the 
applicability of a single (albeit iV-dependent) relaxation 
time model, at least for monodisperse systems. 

C. Coherent strain-activated relaxation 

The next step in constructing a useful microscopic the- 
ory is prediction of t. If one supposes that chains under 
active deformation at strain rate e relax coherently, and 
that the relaxation is strain-activated, r is reduced by a 
factor of Ne relative to its quiescent value. 

We assume relaxation on large scales is coupled to seg- 
mental relaxation. For an a prirori unspecified relaxation 
dynamics in the quiescent state, 

t ~ N^T a (incoherent), (3) 

where r Q is the "alpha" or segmental relaxation time, and 
7 is unspecified and may be TV-dependent. 

Gaussian polymers have chain statistics defined by 
R% = R 2 N, where R 2 is the squared statistical segment 
length. Mathematically, this gives the identity 

dR 2 /dt = NdR 2 Jdt. (4) 

While Eq. @] surely oversimplifies the physics of glasses 
(e. g. it does not hold in quiescent systems because 
R 2 is stationary), nonetheless it suggests an associated 
relaxation rate under active deformation that is N times 
larger (or a time N times smaller) than the value in the 
quiescent state: 

r ^Ary-i Tct (coherent). (5) 

We postulate that Eq. @] becomes valid in actively de- 
formed glasses; segmental rearrangements become "co- 
herent" because rearrangements that restore R 2 towards 
its initial value dominate over those which do not. In 
practice, coherent relaxation is forced by the stiffness of 



the covalent backbone bonds, which have (nearly) con- 
stant length Iq and so maintain (nearly) constant chain 
contour length L = (N — 1)Iq. 

Recent work assists in hypothesizing a more specific 
relaxation dynamics for actively deformed systems. Ref. 
(23| showed that chains in not-too-densely [36| entangled 
model polymer glasses well below T g orient independently 
of one another during active deformation. The behav- 
ior observed was that of individual chains coupled to a 
"mean-field" glassy medium. If relaxations on chain and 
segmental scales are tightly coupled, and chains relax in- 
dependently of one another, 7 = 2 is predicted [2|. Then 

t ~ N 2 r a (incoherent), 

(6) 

t ~ Nr a (coherent). 

Within the framework of Eqs. [31 [SJ and HI when relax- 
ation is strain activated, both r Q [8j and r will be reduced 
by a factor ~ e, so the overall reduction in r during active 
deformation scales as Ne. Note that the scaling analysis 
presented above does not account for additional changes 
in r a arising from other causes, e. g. increased mobility 
associated with yield. Also note that the above argu- 
ments assume e e // and e have the same sign. If deforma- 
tion is reversed, e. g. in a Bauschinger-effect experiment 
[H,|43|, and the sign of e is opposite that of e e //, active 
deformation may not produce coherent relaxation. 

D. Bead-spring simulations 

Some of the arguments made in Sections III BIIII Cl were 
heuristic, and several assumptions were made, so it is 
important to compare the theoretical predictions with 
results from simulations. For example, observations of 
sharp changes in segmental relaxation times with strain 
and significant dynamical heterogeneity 0, IH 0] during 
deformation are seemingly at odds with our postulated 
single, constant relaxation time r. Further, entangled 
chains may not be able to relax coherently if the entan- 
glements concentrate stress, so the reduction of t in ac- 
tively deformed entangled systems may be weaker than 
predicted above. 

The basic ideas presented above can be tested us- 
ing molecular dynamics simulations of the Kremer-Grest 
bead spring model [37| . Polymer chains are formed from 
N monomers of mass m. All monomers interact via the 
truncated and shifted Lennard-Jones potential Ulj — 
Au Q [(a/r) 12 - (0/r) 6 - ((a/r c ) 12 - (a/r c ) 6 )]. Here r c = 
1.5a. Covalently bonded monomers additionally inter- 
act via the FENE potential U FE ne = -(kE$/2) ln(l - 
(r/i?o) 2 ); the canonical (3?| values k — 30u /a 2 and 
i?o = l-5a are employed. All quantities are expressed in 
terms of the intermonomer binding energy uq, monomer 
diameter a, and characteristic time txj = \J ma 2 /u . 
The equilibrium covalent bond length is lo = .96a and 
the Kuhn length in the melt state is Ik — 1.8a. 
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All systems have N c h chains, with N c hN ~ 2.5 x 
10 5 . Periodic boundary conditions are applied along all 
three directions of the simulation cell, which has periods 
L x , L y , L z along the x, y, z directions. Melts arc equili- 
brated 38] and rapidly quenched (k B T = —.Q02u q /t lj ) 
into glasses at T = 0.2u n /ftB ~ 0.6T g . Uniaxial-stress 
compressive deformations are then imposed, using the 
same protocols employed in previous work [22l |: a con- 
stant true strain rate e = L z /L z = — 10~ 5 /tlj is applied, 
with A = L z /L z . A Langevin thermostat with damping 
time IOtlj is used to maintain T, and a Nose-Hoover 
barostat with damping time IOOtlj is used to maintain 
zero pressure along the transverse directions. The val- 
ues of ]e| and T employed here lie within ranges shown 
[39l |40| to reproduce many experimental trends |lj , such 
as logarithmic dependence of a on e and linear scaling of 
the hardening modulus with the flow stress. 

E. Chain conformations under deformation and 
constant-strain relaxation 

Constant strain-rate deformation and constant strain 
relaxation are two of the most commonly performed me- 
chanical experiments. In a constant strain rate exper- 
iment, assuming t is independent of e, i. c. assuming 
polymer glasses are linearly viscoplastic, the solution to 
Eq.dis gU 

£ e//( £ ) = er (1 - exp(-e/er)) = er (1 - exp(-t/r)) . 

(7) 

If a system is deformed to a strain e° and effective strain 
£®ff, and then deformation is ceased, Eq. [5] has the solu- 
tion 

ee//=e° // exp(-t/r). (8) 

In this case, our model predicts slowdown in relaxation 
upon cessation of deformation; the r in Eq. [8] is iVe -1 
times larger than the r in Eq. [JJ 

Figure [5] shows theoretical predictions of Eqs. [JJ5] for 
evolution of e e / / (i) in systems compressively strained to 
e = —1.0 at constant rate (for < t < e _1 ), and then 
allowed to relax at constant strain. Results are plotted 
against t = et, where e is the strain rate applied during 
compression. The solid lines assume 7 = 2, and relax- 
ation is coherent (t cx e^ 1 ^ 7 " 1 ) during active deforma- 
tion and becomes incoherent (r oc AT 7 ) when deformation 
is ceased, consistent with our theoretical picture. For the 
purpose of contrast, the symbols show predictions assum- 
ing that relaxation remains coherent (i. e. r increases by 
only a factor e _1 ) after deformation is ceased. N = 500 
chains orient nearly affinely during strain: e e // — it, 
while short chains orient much less. The solid lines are 
far more consistent with both experiments and sim- 
ulations [f| 0] , which show relaxation slows dramatically 
upon cessation of active deformation. 

Figure [3] shows simulation data for e e f t (t) under the 
same procedure of compression followed by constant- 
strain relaxation. Solid lines for t < 1 are fits to Eq. [7] 




t 



FIG. 2: (Color online) Compressive deformation followed by 
relaxation at constant strain. Curves from top to bottom are 
for N = 500, N = 40, N = 10, N = 5, and N = 3. t = \i\t 
is time scaled by the strain rate applied for < t < 1. For 
t < 1 curves show predictions of Eq. [7] while for t > 1 curves 
and symbols show predictions of Eq. [8] Solid curves assume 
r oc TV 7-1 during deformation and iV 7 at constant strain, 
while symbols assume r oc N J ~ 1 at all t. Here 7 = 2. 

Values for r from these fits are given in Tablelll The data 
are quantitatively consistent with t oc N 1 ^ 1 and 7 = 2. 
In particular, values for t/(N — 1) are nearly constant. 
For this model the entanglement length is N e ~ 85 [43j], 
so this trend spans the range from unentangled to well- 
entangled chains. The effect of entanglements is more 
consistent with an increase in the prefactor of rj (N — 1) 
than a change in 7. 




FIG. 3: (Color online) e e // vs. t for uniaxial compression fol- 
lowed by constant strain relaxation; comparison of theory and 
simulation results. Symbols from bottom to top are bead- 
spring simulation results for iV = 500, 36, 18, and 12. For 
t < 1, lines are fits to Eq. [7J the fit values of r are given in 
Table [TTJ For I > 1 the lines are given by Eq. [8] no further 
fitting is employed, and the values of r from Table [TT] are mul- 
tiplied by iVe -1 , consistent with the transition from coherent 
to incoherent relaxation assumed to occur upon cessation of 
deformation. 

As shown in Table HH the agreement of r with the pre- 
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diction r Q ~ .le _1 under active deformation [8J, where 
we identify t q = t/(N— l) 7-1 , is quantitative. The data 
in Fig.|3]is consistent with our hypothesis that r remains 
constant during deformation; only small deviations from 
the fits to Eq. [7] are apparent on the scale shown. Given 
the large variation in stress as strain increases (see Sec- 
tion III G[) , it is remarkable that a single relaxation rate 
theory fits e e tf so well. However, this is consistent with 
the enhancement of dynamical homogeneity and narrow- 
ing of the relaxation spectrum under flow observed in 
similar models in Refs. 013. In developing our theory 
for the mechanical response, we will assume for simplicity 
that r is indeed independent of e. 

TABLE II: Values of r obtained by fitting e e //(t) to Eq. \7\ 
(with 7 = 2) as a function of chain length N, for flexible bead 
spring chains. Statistical uncertainties on the fits, i. e. the 
error bars for r, increase with TV, to ~ 5% for N = 500. The 
factor of N — 1 in the rightmost column arises because in a 
discrete-bead model, relaxation is associated with the number 
of covalent bonds rather than monomers. The data is from 
Fig. E] Note e = 1Q~ 5 /r LJ so er is of order unity for chains 
with N ~ 10. For this model N e ~ 85. 

N t/tlj t/((N - 1)tlj) 



12 


1.15 


10 5 


1.05 


10 4 


18 


1.78 


10 5 


1.05 


10 4 


36 


4.49 


10 5 


1.09 


10 4 


71 


7.19 


10 5 


1.03 


10 4 


107 


1.12 


10 6 


1.05 


10 4 


250 


2.76 


10 6 


1.11 


10 4 


500 


7.14 


10 6 


1.43 


10 4 



After cessation of deformation, simulation results in 
Fig. [3] suggest Eqs. [5] and [5] give a qualitatively valid 
description of large scale chain relaxation at constant 
strain. The small initial decreases in |e e //| at times 
(i = 1 + 5, i < 1) are larger than the theoretical predic- 
tions of Eq. [5] for small N. This may be attributable to 
chain end or segmental-relaxation effects not included in 
our model. At larger t, however, it is clear that the in- 
crease in r upon cessation of deformation is qualitatively 
captured. 



F. Microscopic theory for viscoplastic stress 

Here we will derive a theory of Neohookean viscoplas- 
ticity for the iV-dependent stress-strain response. The 
work to deform the system is broken into two compo- 
nents corresponding to an isotropic resistance to flow and 
an anisotropic resistance to chain orientation. This gen- 
eral approach has been employed many times before, e. 
g. in Refs. [lj], 53 ■ We present a novel microscopic pic- 
ture that quantitatively associates Gaussian/Neohookean 
strain hardening with the energy dissipated in local seg- 
mental hops, and quantitatively associates the smaller 



hardening observed for shorter chains to strain-activated 
relaxation of large-scale chain orientation (i. e. A ^ A e //). 

From the second law of thermodynamics, the work W 
required to deform a polymer glass is W(e) — AE(e) + 
AQ(e), where E is the internal energy and Q is the por- 
tion of work converted into heat, including both dissi- 
pative and entropic terms. Without loss of generality, 
this can be rewritten as W(e) = W\{e) + IF^e), where 
W2 is the viscoplastic component of the work. W\ cap- 
tures "everything else", such as elastic terms pre-yield, 
and (in principle, though they are not treated in this pa- 
per) other effects such as softening, anelastic energetic 
stresses, and entropic stresses. In experiments, ductile 
deformation of glassy polymers occurs at nearly constant 
volume [1]; thus we assume constant volume deformation 
and treat Wi and W2 as intensive quantities. 

Both rubber elasticity and Neohookean elasticity asso- 
ciate strain hardening with W\. However, experiments 
and simulations [l9|, [22[ have shown that W2 is the dom- 
inant term for strains ranging from the beginning of the 
plastic flow regime to the onset of dramatic "Langevin" 
[17j hardening. The latter occurs at very high strains 
for most synthetic polymers [ll| and has been associated 
with the increase in energy arising from chain stretching 
between entanglements [22|, IHj]. Simulations have pro- 
vided strong evidence that in this regime, W2 is closely 
connected with the same local interchain plastic events 
that control the flow stress |22]. These events have a 
characteristic energy density ~ uo/a 3 , where uq is the 
energy scale of secondary (i. e. noncovalent) interactions 
[49l | . Here we treat the regime where W2 dominates, and 
neglect W%, i. e. we make the approximation W = W2- 

W may be further broken down into "segmental" and 
"polymeric" terms: W — W s 4- W p , where W s accounts 
for the plastic flow stress in the absence of hardening, 
and W p accounts for the viscoplastic component of strain 
hardening. Since this paper focuses on strains well be- 
yond yield, for convenience we choose a standard viscous- 
yield term 

W s = ^(e + e y exp(-e/e y )), (9) 

where e y is the yield strain. 

The stress is given by a — dW/de. It can similarly be 
written as a sum of segmental and 'polymeric' contribu- 
tions 

a = o s + o p = dW s /dt + dW p /de. (10) 

Ref. [22j showed a oc R pi where R p is the rate (per 
unit strain) of plastic events identified by local rearrange- 
ments. The natural correlation length scale for the local 
plastic rearrangements, since according to the above ar- 
guments they are coherent, is R c , and the associated vol- 
ume is V = i? 3 /6 3 / 2 . Here factors of V6 arise from the 
standard relation for the radius of gyration R g of Gaus- 
sian polymers, R g — R c j v6. 

We will now associate polymeric strain hardening with 
W p and replace A by X e ff in order to calculate an A^- 
dependent a p . We postulate that W p is controlled by the 
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increase in V; in other words, strain hardening occurs 
because the volume V controlling W p increases faster 
than it can relax Studies of bidisperse mixtures 

have provided strong evidence that the evolution of A e / / 
(and hence strain hardening) can be understood in terms 
of single chains interacting with a glassy mean field [23[ , 
so it is natural to assume the plastic events are "unary" 
(in the sense that > 2-chain effects are unimportant). 
W p will then scale linearly with p cr , where 



p cr = V6pl /(NR C ) 



(11) 



is the density of coherently relaxing contours, and p is 
monomer number density. Then, incrementally 

AW P ~ (u /a 3 )A{p cr V) = {u /a?)N- l phA{Rl)/&. 

(12) 

Recall R 2 = 3WifTV<?(A e //), where Ik is (here) the "equi- 
librium" Kuhn length in the undeformed glass. Then for 
a deformation increment AA, 



A(R 2 c ) = 3l l K N A(s(A)| 



(13) 



The term in parentheses indicates the difference is eval- 
uated at A = A e //- 

Combining equations [T2TIT31 gives 



AW p = (u /a 3 )pl 2 l K A(g(\ eff ))/2. 



(14) 



This result combined with Eq. [9] gives a prediction for 
the stress: 

°"0e//) = (l - exp(-e/e a ) + pl 2 l K \g(e eff )\/2) , 

(15) 

where g(e) = (3/2)dg/de. The absolute value |<ji(e e //)| 
appears because we have so far treated a as positive. 

Eq. [15] has several interesting features, which we now 
relate to previous models. First, it predicts that a p (at 
fixed strain) increases with increasing Ik- This is con- 
sistent with the well-established result that straighter 
chains are harder to plastically deform 0, [Tl|- However, 
the power of Ik on which a p depends is sensitive to our 
theoretical assumptions, specifically the definition of p cr 
(Eq.[TT]). Other definitions can produce Gr oc l^ 2 or l K , 
but choosing which one is "best" [54| requires greater 
knowledge of the variation in local plasticity (at a micro- 
scopic level) with Ik/Io than is currently available. We 
test Eq. [15] using MD simulations in the following sec- 
tion. While only one value of Ik /a is considered here, 
it is large enough (1.8) that viscoplastic contributions to 
a scaling as l K would be much larger than contributions 
scaling as Ik- 

Second, for long chains, our theory predicts little relax- 
ation during deformation, and a strain hardening modu- 
lus Gr ~ pIqIk- This value is much closer to the effec- 
tive constraint density measured in (NMR) experiments 
of deformed glassy samples [14] than to the entropic pre- 
diction Gr = p e ksT. It must be noted that increasing 
Ik also increases p e 52]. However, considerable evidence 



(e. g. d, [H [H, [H, S3, [H HU) suggests that chain ori- 
entation and local, secondary interactions which act over 
scales ~ (Iq ~ a ~ Ik) are the true controlling factors for 
<7, at least during the initial stages of hardening. Finally, 
Haward postulated that Gr arises from the contraints 
imposed by the mesh of uncrossable chains fill ]; our ar- 
gument that hardening scales with consistent with 
this hypothesis. 



G. Predictions for stress-strain curves 

In Equation [15] we have assumed that all contribu- 
tions to a scale with a single energy density (i. e. stress) 
(uo/a 3 ). This is an approximation, since flow and hard- 
ening stresses have been observed to be linearly rather 
than directly proportional 4l|. However, the 'constant 
offset' term in this linear relationship is often fairly small 
compared to the linear term (especially for T well below 
T g [H, Hi)]), so the approximation is reasonable. There- 

and scale it out. This 



fore we associate uo/a 3 with cr// 
gives 



a V e eff) = a — — . 



(16) 



where A is a prefactor arising from our neglect of pref- 
actors in the above analysis. Note that all "thermal" 
aspects of our theory are implicitly wrapped into A and 
r a . In both real and simulated systems a tin, ,, and A are 
approximately proportional to (1-T/T g ) tUM IM 113 . 
so this scaling should remove much of the T-dependence. 
The merits of "multiplicative" forms like Eq. rTB] for 
predicting stress have been discussed recently in Refs. 
[43, More sophisticated models (e. g. Ref. [§]) ex- 
plicitly treat the variation of A with e and T, and/or the 
variation of r a with a, T, and local microstructure. Here, 
however, A and (for fixed N) r Q are treated as numerical 
constants. 

To predict the TV-dependence of u(e), an analytic form 
for g(e e ff(e)) is obtained by plugging in the solution for 
e e //(e) from Eq. [7] into the form of g(e) for uniaxial de- 
formation (Table [H with exp(e) = A): 



S(e e //(e)) = exp[2er(l - exp(-e/er))] 
— exp [— er(l — exp(— e/er))] . 



(17) 



Figure [4] shows predictions of Eqs. IT&HTTl for er*(e) in uni- 
axial compression at various N, with A = 1, \e y = .02], 
p = 1.0a -3 , Iq = .96a, and Ik = 1.8a; the latter three 
are chosen to match the flexible bead spring model em- 
ployed in the simulations. Solid lines assume r oc N J ~ X 
with 7 = 2 as discussed above, while symbols assume 
coherent chain relaxation is not important and r cx N 1 . 
The solid lines are qualitatively consistent with simula- 
tions [13, H3, while the symbols are inconsistent. Both 
show increasing strain hardening with increasing TV, but 
in the latter case hardening increases much faster and 
saturates at a much lower value of N than is realistic. 
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FIG. 4: (Color online) Stress-strain curves predicted by Eqs. 
1161171 with A — 1. —a* and — e are shown because stress 
and strain are negative for compression. Curves from top 
to bottom are for TV = 500, N = 10, N = 5, and TV = 3. 
Predictions in the TV — > oo limit are not distinguishable from 
TV = 500 predictions on the scale of this plot. Solid curves 
assume er = 73 TV 7 " 1 , while symbols assume er = 73 TV 7 . In 
both cases, 7 = 2, and as suggested by Ref. [g], B — 0.1. 



For example, the r oc N 1 predictions for TV = 40 and 
TV = 500 are indistinguishable on the scale of the plot. 
Results for tension are not presented here because our 
model includes no asymmetry between tension and com- 
pression 



58|. 



The large-strain (e 3> t y ) mechanical response pre- 
dicted by our model varies continuously from perfect- 
plastic flow (cr*(e) — > a constant (jfi ow ) to network- like 
polymeric response (c*(e) — > 1 + plglKg(^}) as er varies 
from zero to oo (equivalently, as TV increases). In be- 
tweeen these limits, the response is "polymeric viscoplas- 
ticity". Note that since we treat e e // as microreversible, 
a more accurate term for our model is "viscoelastoplas- 
ticity" [1^]. However, since our model does not treat 
stress relaxation after cessation of deformation, a* may 
be regarded as a (reduced) orientation-dependent plastic 
flow stress. 

We now compare theoretical predictions for stress- 
strain curves to results from bead-spring simulations. In 
FigureOH dashed lines show predictions of Eqs. lTolfTTl with 
A = 0.43, values of r taken from Table [TTl and the same 
values of |, l and Ik as in Fig. 3J The value of A ob- 
tained in Fig.[5]is comparable to (1 — T/T g ). Stress-strain 
curves from simulations are shown as solid lines. Note 
that stresses and strains are negative in compression. 
Panel (a) shows a*, while panel (b) shows its dissipative 
component <r Q * = (a* - (u / a 3 )' 1 dU / de) [||. While 
this definition of a®* includes entropic terms, these are 
known to be only of order 1% of the total stress at this 
T [H. 

In both panels, the correct trends are predicted, and 
quantitative agreement is within ~ 20%. Very short 
chains (TV = 4) show nearly perfect-plastic flow, while 
longer chains show strain hardening, similar to results an- 
alyzed in many previous studies. The quantitative difler- 
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FIG. 5: (Color online) Stress-strain curves (a*) from bead 
spring simulations (solid lines) and theoretical predictions of 
Eqs. DHTI] [(a) cr*, (b) ° Q *]. Values of TV from top to bottom 
are 500, 36, 18, 12, and 4. The top 4 are from the same 
simulations shown in Fig. [3] and use the fit values of r (Table 
HI)), while the bottom theory curve assumes er = .1(TV — 1). 



ences at small strains arise primarly from our oversimpli- 
fied treatment of yield At large strains, comparing 
panels (a-b), it is apparent that differences between pre- 
dictions of Eqs. fTBlfTTl and simulation results arise largely 
from energetic terms associated with strain hardening, 
i. e. covalent bond energy and additional plastic defor- 
mation arising from chain stretching between entangle- 
ments. Experiments on amorphous polymer glasses (e. g. 
[ll|) have shown a similar percentage of the total stress is 
associated with energetic terms, strain softening, etc., so 
the agreement between theoretical predictions and bead- 
spring results is satisfactory given the simplicity of our 
model. In panel (b), the "spike" in simulation results 
at small strains for V > 4 reflects yield and subsequent 
strain softening. 



H. Nonaffine displacement and plastic deformation 

We have argued above that coherent relaxation is 
driven by the stiff covalent bonds and the need to main- 
tain constant chain contour length L. Since long chains 
have er ^ 1 and deform affinely on the end-to-end scale, 
they must deform nonaffinely on smaller scales to main- 
tain chain connectivity. It is interesting to relate this 
nonaffine deformation to the chain stretching that would 
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occur if deformation were affine on all scales and L was 
not constrained, i. e. the limit of a Gaussian coil with zero 
spring constant 0] embedded in a deforming medium. 
For uniaxial tension or compression the nonaffine dis- 
placement should be given by D^ a ~ 5(A) — 1, with 

-w-SfHH^)- (18) 

Figure [5] shows data for the squared nonaffine displace- 
ment of monomers, D^ a = {(r — Ar*o) 2 ), where the 
monomer positions are {r} at stretch A and {rb} in the 
initial state, for N = 500. Data for T = .Oluo/^s is 
shown to minimize the thermal contribution to D„ a . A 
fit to 5(A) — 1 is also displayed. There is qualitative 
agreement at large strains (A -C 1), and the underesti- 
mation of Df la at smaller strains is largely attributable 
to smaller scale (i. e. incoherent) plasticity on scales ap- 
proaching the monomer diameter a. Another reasonable 
form for fitting to D\ a is g(X), which would suggest D\ a 
scales with W p at large strains. This form gives slightly 
less good fits to our data, but in any case the principle 
shown in Fig. [5] is the same as outlined in Ref. (6(j; for 
long chains, nearly affine deformation at large scales (in 
our language, A ~ A e //) increasingly drives nonaffine dis- 
placements (i. e. plastic activity) at smaller scales, lead- 
ing to strain hardening. This effect weakens and values 
for both a and D\ a decrease (for large strains) with de- 
creasing N. 
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FIG. 6: (Color online) Nonaffine displacement at low T in 
systems of long chains. Circles show data from bead-spring 
simulations (with TV = 500 and T = .01«o/fcs). The solid line 
shows a fit to D 2 na / a 2 = C(S(X) - 1) (see Eq. [18]) with C = 
7.39. Note that the maximum value of Dn„ is well below the 
squared tube diameter d^/a 2 ~ 100 @, HJ, consistent with 
the picture that entanglements do not control the mechanical 
response of these systems. 



III. DISCUSSION AND CONCLUSIONS 

We derived a simple theory for polymeric strain hard- 
ening based on the notion that the increase in stress be- 



yond yield amounts to an increase in the "flow" stress 
in an increasingly anisotropic viscoplastic medium. In 
the strain hardening regime, long and short chains relax 
on large scales via similar mechanisms (r cx e N) and 
longer chains show greater hardening because they can- 
not relax on large length scales over the timescale of the 
deformation. In practice, this is due to the high inter- 
chain friction in the glassy state. 

We theoretically predicted, and provided evidence us- 
ing simulations, that coherent chain relaxation driven by 
resistance to chain contour length increase is a key factor 
in the large-strain mechanical response of uncrosslinked 
polymer glasses. Coherent relaxation reduces the domi- 
nant (chain scale) relaxation time r by a factor Ne during 
active deformation. Further, we claim that the increase 
in relaxation times when deformation is ceased is at least 
partially associated with the fact that relaxation need no 
longer be coherent. Our results are consistent with many 
previous simulations and experiments, and semiquanti- 
tatively capture the increase in strain hardening as chain 
length increases. 

Much recent work has emphasized the predominantly 
viscous/viscoelastic nature of polymeric strain harden- 
ing, at least prior to entanglement-stretching. Recent 
experiments and modeling [471 . l53l . 1571 l64j have provided 
strong support to the notion that entanglements play 
only a secondary role in glassy-polymeric strain harden- 
ing, at least for the majority of synthetic polymers and 
in the weak hardening regime. Refs. [HI, [64[ also empha- 
sized the role of mcltlikc relaxation mechanisms in the 
marginally glassy state. 

The present work, which focuses on relaxation mech- 
anisms, is consistent with this these trends. We quan- 
titatively related strain hardening to local plasticity at 
the segmental scale, and showed that the power-law de- 
pendence of the large-scale chain relaxation time t on N 
during active deformation, r cx A 7-1 , is consistent with 
7 = 2, the same value as the Rouse model for unentan- 
gled polymer melts @, HH . Another "Rouselike" aspect 
is the apparent tight coupling of r to a single microscopic 
relaxation time r a . 

X e ff can now be accurately measured in scanning 
near-field optical microscopy (SNOFM) experiments [65j . 
which have shown that e e f t < e for entangled chains de- 
formed slightly above T g . Analogous studies, well be- 
low T g , could be performed to test the theory devel- 
oped here; modern neutron scattering techniques might 
be employed for the same purpose. Additionally, defor- 
mation calorimetry (DC) experiments [66j can be per- 
formed to better understand the dissipative contribution 
to the stress. We are not aware of any studies in which 
the SNOFM or DC methods have been applied to study 
strain hardening in polymer glasses. 

Our model is minimal. It cannot quantitatively predict 
stress-strain curves either at small or at very large strains, 
because it neglects energetic components of stress and 
also strain softening, which play important roles in glassy 
polymer mechanics. However, these limitations do not 
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violate the spirit of our modeling effort, which was to 
illustrate the role of coherent relaxation in controlling 
large-scale chain conformations and influencing stress in 
actively deforming polymer glasses. 

The theory presented here serves as a complement to 
a recent microscopic theory by Chen and Schweizer 
which also provides a unified description of plastic flow 
and strain hardening in polymer glasses. Ref. |9j assumes 
X e ff = A, and thus appears to neglect an important re- 
laxation mechanism. Also, Ref. [9j is based on an ex- 
tension of liquid state theory to the glassy state [|| and 
(formally) breaks down at zero temperature. In con- 
trast, the theory developed here assumes relaxation is 
dominated by strain activated processes, and so should 
be most accurate at low temperatures; we expect it to 
break down concurrently with the validity of the "mean 
field" , independent-chain-relaxation behavior [23[ of A e / / 
as T -> T g . 

An ideal approach would account for both thermally 
and strain activated relaxation without sacrificing sim- 
plicity. For example, r should have explicit temperature 
dependence, and thermal activation should produce log- 
arithmic corrections to e _1 scaling. The dependence of r 
on \<j\ (from stress-assisted rearrangements) should also 
be captured. Combining our theory with a more sophis- 



ticated theory for a fi ow 

(e. g. Refs. Q or 0) might be 
fruitful. Alternatively, modern techniques in nonequilib- 
rium thermodynamics may prove a useful tool for treat- 
ing finite-T effects. 

Clearly, further work is necessary to quantitatively pre- 
dict stress-strain and stress-relaxation curves. In par- 
ticular, improved understanding of the effects of chain 
stiffness is required. It seems certain that microscopic 
structural detail at the Kuhn scale (e. g. chemistry- 
dependent effects) exerts significant influence on segmen- 
tal relaxation processes, and it is probable that these 
effects couple to chain-scale relaxation (see e. g. Refs. 
[1, HH, Hl|, |60|, [fi3)- Studies using more chemically realis- 
tic models, or real polymers, would be welcome. 

All MD simulations were performed using LAMMPS 
[68| . Mark O. Robbins contributed significantly to this 
project, with numerous helpful discussions. Kenneth 
S. Schweizer, Grigori Medvedev, Kang Chen, Daniel J. 
Read, and Edward J. Kramer also provided helpful dis- 
cussions, and K. S. provided the original concept for \ e ff- 
Gary S. Crest provided equilibrated TV = 500 states. 
The authors wish to thank the KITP Glasses '10 confer- 
ence, which was the setting of many of these discussions. 
Support from NSF Awards No. DMR-0520415 (RH) and 
DMR-0835742 (RH, CO) is gratefully acknowledged. 
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